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Abstract 

In spite of its undeniable success, there are still open questions regarding Tsallis' non-extensive 
statistical formalism, whose founding stone was laid in 1988 in JSTAT. Some of them are concerned 
with the so-called normalization problem of just how to evaluate expectation values. The Jaynes' 
MaxEnt approach for deriving statistical mechanics is based on the adoption of (1) a specific en- 
tropic functional form S and (2) physically appropriate constraints. The literature on non-extensive 
thermostatistics has considered, in its historical evolution, four possible choices for the evaluation 
of expectation values: (i) 1988 Tsallis-original (TO), (ii) Curado-Tsallis (CT), (iii) Tsallis-Mendes- 
Plastino (TMP), and (iv) the same as (iii), but using centered operators as constraints (OLM). 
The 1988 was promptly abandoned and replaced, mostly with versions ii) and iii). We will here 
(a) show that the 1988 is as good as any of the others, (b) demonstrate that the four cases can 
be easily derived from just one (any) of them, i.e., the probability distribution function in each of 
these four instances may be evaluated with a unique formula, and (c) numerically analyze some 
consequences that emerge from these four choices. 

PACS: 05.30.-d, 05.30.Jp 
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I. INTRODUCTION 



James Bernoulli published in the Ars Conjectandi the first formal attempt to deal with 
probabilities in 1713 and Laplace further formalized things in his Theorie analytique des 
Probabilites of 1820. As we know, in the intervening years Probability Theory (PT) has 
grown into a rich, powerful, and extremely useful branch of Mathematics. Contemporary 
°hvsics heavily relies on PT for a large part of its basic structure, Statistical Mechanics 
[J of course, being a most conspicuous example. One of PT basic definitions is that 
of the mean value of an observable A (a measurable quantity), which in its original form it 
is taken as a lineal combination of the different probabilities that A adopts for the different 
accessible states (we denote them by the index i). Let A stand for the linear operator or 
dynamical variable associated with A. Then 

(A) = Y,P*A. (1) 

i 

It is well known that, in some specific cases, it becomes necessary to use "weighted" mean 
values, of the form 

(A) = /(ft) A, (2) 

i 

with / an analytical function of the pi. This happens, for instance, when there is a set of 
states characterized by a distribution with a recognizable maximum and a large tail that 
contains low but finite probabilities. One faces then the need of making a pragmatical 
(usually of experimental origin) decision regarding /. 

Nonextensive thermostatistics (NET) 



Tsallis' nonextensive thermostatistics (NET) is nowadays considered by many authors to be 
a new paradigm for Statistical Mechanics (SM) Q, ff. It encompasses a whole family of 
SM-formulations, all of them based upon Tsallis' nonextensive information measure 



I Ml 

4, 



1 _ v q 

S q = k B l ^=i p ' , (3) 
q — 1 

where ks stands for Boltzmann's constant and {p n } is a set of normalized probabilities. The 
real parameter q is called the index of nonextensivity. The literature on Tsallis' thermo- 
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statistics considers three possible choices for the evaluation of mean- or expectation-values 
(EVs) within the nonextensive scenario. Since EVs are always regarded as constraints in the 

n 

associated MaxEnt approach [6j, three types of q— statistics will thereby ensue. Consider 
the physical quantity O that, in the microstate n (n = 1,2, ...,W), adopts the value o n . 
Let p n stands for the microscopic probability for such microstate. The O— EV is evaluated 
according to three distinct recipes, that we illustrate below with reference to Gibbs' canon- 
ical ensemble, where it is assumed that the a priori physical information is given by the 
mean value of a Hamiltonian (O = H) of spectrum {oj = ej}, and there are two Lagrange 
multipliers, one associated to normalization (a) and the other to the mean energy (inverse 
temperature ft). 

i) First choice: Tsallis (1988) Q 

w 

(0) = J2p^- (4) 

i=l 

The associated, Tsallis canonical distribution was (see the pertinent demonstration below 
with reference to Eq. (|21j) ) 

W = [{l-qXa + fteJ/q} 1 ^. (5) 

The awkward location of the multiplier a makes Pi hard to evaluate, save for very special 
cases like that of the harmonic oscillator [7(. This motivated people to look for better 
alternatives (see below). Other normalization choices thus emerged that typically employ 
"weighted" mean values and avoid the explicitpresence of the multiplier a: 

ii) Second choice: Curado - Tsallis (1991) 8] 

w 



(0> = 5>fr' (6) 



i=i 



with a canonical distribution given by 



p t = (Z,)- 1 [l-(l-q)ft 6i ] 1 ^ 

Z q = J2ll-(l-<l)Pzi} m - q) - (7) 

i 

iii) Third choice: Tsallis - Mendes - Plastino (1998) 



We need to define here a normalization factor (not a Lagrange multiplier) 



X n 



Erf- 



(9) 



The pertinent canonical distribution, that is self-referential because the pi depends upon X q . 
Since it reads [l| 



Pi 


= z- 1 [1 


X q 









X n 



(l-g)ff 

X n 



Si - U„ 



1/(1-9) 



(10) 



this means that even if one knows U q (as a real number) and the spectrum (we call this the 
canonical practitioner's input), this is not enough to obtain straightforwardly the probability 
distribution. One still needs to repeatedly iterate (JTUj) until numerical self-consistency is 
achieved, which can be difficult in practice. In order to overcome such a problem still 
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which avoids the self- 



another choice can be made, the so-called OLM one 
reference problem. The corresponding probability distribution arises from solving anyone of 
two different variational problems 13], namely, those arising after 



1. replacing S q by Renyi's information measure S R 



R _ ln Li=lP. 



s 



(i 



(ii) 



and keeping the TMP constraints. We call this procedure the TMP-R one. 



2. Maintain Tsallis' entropy, but replace the TMP constraints by "centered" mean values. 
We call this procedure the OLM-one. 



In both cases one finds 



141 



Pl = Z~ l [l - (l-q)(3(e t - U q )} 1 /^ 
X q = (caracteristic identity) 

Z q = J2ll-(l-q)f3(e t -U q )] 1/ ^\ (12) 

i 

and it should be noted that, with access to the canonical practitioner's input, no numerical 
self-consistency treatment is now required. 



Reasons for using weighted mean values 



In the first stage of NET development, its pioneer practitioners made the pragmatic choice of 
using "weighted" mean values, of rather unfamiliar appearance for many physicists. Why? 
The reasons were of theoretical origin. It was at the time believed that, using the familiar 
linear, unbiased mean values, one was unable to get rid of the Lagrange multiplier associated 
to probability-normalization (see below), usually denoted with the letter a in the Litera- 
ture. Since the Tsallis formalism yields, in the limit q — > 1, the orthodox Jaynes- Shannon 
treatment, the natural choice was to construct weighted EVs using the index q 



the so-ca 



i 

led Curado-Tsallis unbiased mean values (MV) 



(13) 



As shown by Plastino and 



Plastino 3], employing the Curado-Tsallis (CT) mean values allowed one to obtain an an- 
alytical expression for the partition function out of the concomitant MaxEnt process jl^ . 
This EV choice leads to a non extensive formalism endowed with interesting features: i) the 
above mentioned property of its partition function Z, ii) a numerical treatment that is rela- 
tively simple, and iii) proper results in the limit q —>■ 1. It has, unfortunately, the drawback 
of exhibiting un-normalized mean values, i.e., ((l)) q ^ 1- The latter problem was circum- 
vented in the subsequent work of Tsallis-Mendes-Plastino (TMP) 
CT treatment by employing mean values of the form 



9[, that "normalized" the 



(14) 



Today, this is the orthodox normalization procedure for NET practitioners. However, the 
concomitant treatment is not at all simple. Numerical complications often ensue, which has 
encouraged the development of different, alternative lines of NET research. See, for instance 



(a by no means exhaustive list), |l6l. Il7 



19 



20 



21 



and references therein. 



Present goals 



In this work we revisit the MaxEnt approach with unbiased mean values as prior knowledge, 
and show that it is possible to obtain a probability distribution that i) does not explicitly 
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contain the obnoxious a factor, ii) allows a seemingly analytical expression for the partition 
function out of the concomitant MaxEnt process (although it remains self-referential), iii) 
it is normalized, and iv) leads to the proper result in the q — > 1 limit. It does exhibit 
some numerical drawbacks, similar to those faced by the TMP treatment, but without the 
theoretical limitations that affect both the CT and the TMP formalisms. Since the treatment 
here advanced allows one to avoid weighted expectation values, it might be preferable to 
both the CT and and TMP ones on Ockham's razor-grounds. 

Moreover, we will show that all four choices described above, namely, (jU, (JHJ), ©, and 
OLM, are equivalent in the sense that the four cases can be easily derived from just one 
(any) of them, i.e., the probability distribution function in each of these four instances 
may be evaluated with a unique formula, that (in the canonical instance) reads, for a given 
q— value 

[l-(l-g*)^] 1/(1 - g * } 

[l-a-g*)^ 117 ^ 

for (3*Ei < — , (15a) 



Pi 



exp g {-(3* Sj) r ^ 1 

Pl = 0; if (3*6, > — L;, (15b) 
1 — q* 

where f3* is an appropriate pseudo-inverse temperature (= l/ksT* ) which depends on (3, 
exp g is the well-known q— generalized exponential ^| 

exp ? (x) = [l + (l-g)x] 1/(1 " 9) , (16) 

and q* is an effective nonextensivity index, suitably related to q (see Table 1). The extension 
of (|15|) to the general case of M observables Oj , (j = 1 , . . . , M) , that involves M Lagrange 
multipliers Xj, (j = 1, . . . , M), is straightforward. 



II. RESCUING TSALLIS 1988 VERSION 
Classical canonical treatment 



For didactic purposes we will restrict ourselves first to the canonical ensemble and the 
classical formulation. The forthcoming section will treat things in quantal fashion for general 
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ensembles. The Boltzmann constant will be taken as unity in this Section. We write the 
internal energy (the mean value of the Hamiltonian) as the linear probabilities combination 



w 



(17) 



i=i 



where the sum is taken over the W micro-states of the system. The MaxEnt procedure entails 
maximizing Tsallis information measure subject to the constraints given by the internal 
energy and normalization of the probabilities 



5> = L 



The quantity to be maximized is 



(18) 



L T = S q - a^pi - pU g , 



i.e., 



dL 



T 



Q q-l 



dpi 



Pi 



a- f3Hi = 0, 



wherefrom we easily find 



Pi 



q) 



q 



a + f3Hi 



1/(9-1) 



(19) 



(20) 



(21) 



This is the result obtained and used in the very first paper on the Tsallis formalism [7], and 
subsequently abandoned i) because of the presence of the a term in (}2"T]) and ii) due to the 
fact that it was seemingly impossible to obtain a proper partition function. The step that 
transforms the above treatment into a useful one is performed via a simple trick. Note that 
we can obtain a by first multiplying Eq. (J2(Jj) by Pi and then summing over i, which yields 



a 



q 



q 



-X, - (3U q . 



Once we have a, the problem is solved. If we place a into Eq. (|21j) we obtain 



Pi 



q-l 



V(9-l) 



(22) 



(23) 
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where 5 Hi = Hi — U q . To explicitly obtain the partition function is now a simple task. Using 
Eq. (fTSj) one immediately finds 



1 - *±4-8Hi 

q X q 



1 1/(9-1) 



Xc 



1/(1-9) 



so that, if we agree to call 



(24) 



i _ £zl A,ffi 



we obtain the well-known, typical TMP result 



q X q 

a 



1/(9-1) 



(25) 



*9 = ^ 



(26) 



and 



Pi 



1 V(9-l) 



(27) 



It is clear that the need of having a real probability forces one to retain only the positive 
values of the expression between brackets (cut-off condition). Notice that the result (j2HJ) 
resembles the TMP one and that the presence of the quantity X q makes the probability 
distribution a self referential one, although the proper limit when g — f 1 is indeed obtained 

e f3U e -f3H t e -pHi 
^9=1 = e /3U^ e -/3Hi = T re -/3^ • ( 28 ) 

It remains as yet to ascertain whether the entropy is indeed a maximum, a point that has 
not been sufficiently discussed thus far. Taking the second derivative of L? we easily obtain 



d 2 L 7 



9-2 

-QPi , 



(29) 



i.e., for q > the entropy will be a maximum. This result entails that negative q— values are 
not acceptable. 



General quantal treatment of the Tsallis' 1988 version 

Within the general quantum nonextensive framework, all information concerning the system 
of interest is provided by the statistical operator p. To find it one has to extremize the 
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information measure 

S q = (30) 

where now 

X q = Tr , (31) 
subject to the normalization requirement 

Tr(p) = l (32) 

and the assumed a priori knowledge of the expectation values (EV)s of, say, M relevant 
observables Oj with j = 1...M 

(6 1 ) = Tr(p6 3 ). (33) 

The density or statistical operator p is here the quantum version of the set of probabilities of 
the preceding Subsection. Following exactly the same procedure as before, with M Lagrange 
multipliers Xj, we obtain a density operator of the form 

P = (34) 

with a partition function given by 

Z q = Tr /Vfe- 1 ), (35) 

where 

/ = 1 - (q ~ 1) E 7?*°' ( 36 ) 

is sometimes called the configurational characteristics. It is finite only if its argument is 
positive. / is null if its argument is negative (cut-off condition). The quantities 

SOj = 6 3 - {0 3 ) q (37) 

are usually called centered operators. 



Connections between de 1988 and the TMP formalism 

Comparing Eq. (pMjl with the density operator of TMP treatment , we see that it is easy 
to pass from one to other via an appropriate mapping. If the subindex u stands for the 1988 
NET version we have 
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Q-Mu ^ 1 ~Q\tmp 

V<?L ^ Htmp- ( 38 ) 



An interesting consequence of Eq. (}3q) ensues: as S q will be convex for q\ u > we discover 
that the convex region for S q within the TMP formalism will be q\rMP < 2. Now, it has 
been shown in Ref. jlfll. \v\ that in the TMP instance on also finds the requirement q > 0. 
Thus, TMP and 1988 results should be valid for q e (0, 2). 



III. EQUIVALENCE OF THE DIFFERENT TSALLIS VERSIONS 

As stated above one takes, for a numerical treatment of the canonical ensemble i) the 
internal energy U q as a numerical datum and ii) the energy spectrum. This provides one 
with a probability distribution. We have called i) and ii) the canonical practitioner's input. 
We advance here an alternative route (we call it the parametric approach), in which, instead 
of giving U q as a numerical input data, one a priori fixes the numerical value of an effective 
inverse temperature (3*. With such an alternative procedure one is able to achieve two 
different goals, for a given energy- spectrum: 

• a unification of all Tsallis' techniques, from a theoretical viewpoint, and 

• a very convenient numerical procedure, from a practical perspective. 

Indeed, by recourse to an analysis similar to the one given above, effected now for all 4 types 
of Tsallis distribution mentioned in the Introduction, it is easy to show that this parametric 
technique can i) link the four NET versions and ii) give the pertinent four distribution-forms 
a "universal" (common) appearance. The ensuing results are summarized in Table 1 for the 
canonical ensemble. The energy spectrum {Ei} constitutes the basic a priori data. A power- 
law distribution- form (common to all 4 Tsallis versions) of exponent q*(q) can be built up, 
for any q— value, that reads 

Pi = Z~ x exp g */ 9 )(— j3*£i)\ whose (version — dependent) ingredients are : (39) 

• an effective inverse temperature /?*, given a priori as a number, 
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Z q , the normalization constant (partition function), and 



• an "effective" nonextensivity index q* = q*(q), given in the Table. Indeed, save for 
the case of the 1988 version, we have q* = q. 

In all four cases, U q is evaluated a posteriori, as indicated in the Table. One can demonstrate 
that the appropriate entropic form S q is (a posteriori) maximized by (|39|) (for this particular 
Ug). Also, 

1. T is the "true" thermodynamic temperature that a thermometer would measure, and 

2. (3 is the Lagrange multiplier of the putative MaxEnt process one could have employed 
to reach the "parametric" PD (|39|) using the canonical practitioner's input. 



By TO we mean Tsallis' 1988 original treatment, and by TMP-R the OLM one HJ, |ll 
12j. From the canonical practitioner's viewpoint (3* is, save for the CT and OLM (TMP- 
R) instances, a rather formidable quantity that makes the probability distribution a self- 
referential one (the signature of the self-referential character is the appearance of X q in the 
expression for (3*). This is not so from the parametric perspective. As explained below, if 
one fixes (3* at the outset, considering it as an input parameter, the self-referential difficulty 
vanishes. One is then able to numerically build in fast fashion the PD ()39|) for very many 
(3*— values in the interval [0, oo] and then "read" off the Table the appropriate values for 
T(f3), U g , Xq associated to any given (3*. 



Table 1: Non-extensive Statistics 


Probability Distribution Function 






Pi = Z 1 expg, ( 


-(3*8,), i£/3* £i < j^; Z = E^exp^ (-/?*£,) 






Pi = o, 


if/?*£ l>T ^ 




(3* 


q* 


u q 


T 


TO 


(3* 


_ P 


q* = 2 - q 


U T q° = Y. % P^ 




<?X 9 + (</-l)/9C/ ru 


B IXq 


CT 


(3* = 13 


q* = q 


Uq T = J2i Pl £ i 


k B T = k B T* 


TMP 


(3* 


_ P 


q* = q 


JjTMP TriPi^i 


k B T = k ^-^ MP 


Xq + {l-q)f3ljrMr 


i EiPi 


Xq 


TMP-R 


13* -- 




q* = q 


ttTMP-r _ Y, l p q i £i 


k B T = k B T* -(1-q) UT Mp - R 


i+(i^)Pu'r ip - H 


q EiPi 
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Notice that, if any one of the four Tsallis' probability distributions (PDs) is available, we 
can with it, using the appropriate equations of the version at hand, evaluate the quantities 
U q , X q , and from them obtain T. Now for translating into another Tsallis' version we 
-use our PD to evaluate the U q — value of the new version. With it and 
-the above values for X q and T, compute the new (3*, appropriate for this second version. 
- Finally, employ Eq. (ffilj). 

In order to gain some insight onto the properties of each Tsallis- version we have computed, 
for two different q values, namely, q = 0.6 (Figs. 1 and 3) and q = 1.5 (Figs. 2 and 4) the 
relative internal energy XJ q j{hv) (Figs. 1 and 2) and the specific heat C q /kB (Figs. 3 and 
4) for a system of equally spaced levels labelled by a "quantum number" i, like in the case 
of the harmonic oscillator, for instance (inter-level distance = Eq), with i = 0,1,2,... and 
Eq = hp. For all four Figs., the horizontal axis represents k&T '/(hu) and the q = 1 results are 
also shown for comparison. A glance to Figs. 1-4 shows that the four Tsallis versions yield 
different results, and that the OLM instance is the one that resembles the orthodox q = 1 
situation in closer fashion. In particular, notice in Fig. 2 that the OLM— q = 1.5 curve is 
identical to that for q = 1. Specific heat undulations are a consequence of Tsallis' cut-off, 



as explained in Ref. 



12|. The asymptotic behavior when T — > oo for U q and C q is given by 



the following expressions: 



U q ~T° 



C q ~ 



T 



a-1 



where the appropriate a-exponent values, for each Tsallis version, are shown below: 



TO 


CT 


TMP 


TMP -R 


a= l - 

i 


a = 2 — q 


a= l - 

i 


a = 1 



Advantages of the parametric form (j39|) 

It should be remarked that using the probability distribution (PD) ()39|) exhibits the advan- 
tage of avoiding self-reference no matter which of the normalization choices one uses. In 
writing down the PD in the fashion 



Pi 



Z 1 exp ,{-f3*Ei)- } Z q 



E 



exp. 



(40) 
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you can numerically input, for a given 

q* = 2 — q; for the 1988 version 

q* = q; for all the other versions, (41) 

any f3*— value whatsoever and you get the PD. No self-reference here. Since the data are the 
spectrum, (3*, and q*, we speak of a "parametric" calculation. With such a PD you evaluate 
both X and U q and then proceed to ascertain the value of the physical temperature T using 
the relations given in Table 1. From a numerical viewpoint therefore, (J4(Jj) provides one with 
a convenient tool to work out nonextensive thermostatistics problems. 

CONCLUSIONS 

We have revisited the four different extant versions of Tsallis' nonextensive thermostatistics 
and shown that they can be unified (Cf. Table 1) in the sense of deriving the transformations 
that provide a link between any given pair of versions. It must be emphasized, however, 
that numerical results are different for the distinct versions. In particular, we have shown 
that the currently abandoned original Tsallis version of 1988 is a quite legitimate one. 
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FIGURE CAPTIONS 

Fig. 1: Relative internal energy U q /hu vs. k B T/hu (q = 0.6) for the four different 
instances of Tsallis-normalization (see text). 

Fig. 2: Relative internal energy U q /hu vs. ksT/hu (q = 1.5) for the four different 
instances of Tsallis-normalization (see text). 

Fig. 3: Specific heat Cg/ks vs. ksT/hu (q = 0.6) for the four different instances of 
Tsallis-normalization (see text). 

Fig. 4: Specific heat C q /k B vs. k B T/hu (q = 1.5) for the four different instances of 
Tsallis-normalization (see text). 
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